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Abstract 

For sufficiently high Mach numbers, small disturbances on a supersonic vortex sheet 
are known to grow in amplitude because of slow nonlinear wave steepening. Under the same 
external conditions, linear theory predicts slow growth of long-wave disturbances to a thin 
supersonic shear layer. An asymptotic formulation is given here which adds nonzero shear- 
layer thickness to the weakly nonlinear formulation for a vortex sheet. Spatial evolution is 
considered, for a spatially periodic disturbance having amplitude of the same order, in 
Reynolds number, as the shear-layer thickness. A quasi-equilibrium inviscid nonlinear 
critical layer is found, with effects of diffusion and slow growth appearing through a 
nonsecularity condition. Other limiting cases are also considered, in an attempt to determine 
a relationship between the vortex-sheet limit and the long-wave limit for a thin shear layer; 
there appear to be three special limits, corresponding to disturbances of different amplitudes 
at different locations along the shear layer. 

1. Introduction 

According to the linear theory, a vortex sheet in supersonic flow is neutrally stable if 
the Mach number is large enough. For example, if the temperature is the same on both sides 
of the sheet, neutral stability is predicted when the difference in Mach numbers on the two 
sides is greater than 2^2; otherwise disturbances grow exponentially with time. Thus if the 
velocities are equal and opposite above and below, the Mach number in each part of the flow 
must be greater than V2 to avoid instability. For higher Mach numbers, there is a stationary 
mode in this coordinate system, as well as a pair of propagating modes that move at speeds 
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depending on the Mach number. The definitive presentation of these results was given by 
Miles [1], who gives additional references. 

Artola and Majda [2] have considered the nonlinear growth of a disturbance to a 
vortex sheet that is caused by an incident sound wave. For incidence angles such that the 
component of wave number along the sheet equals the wave number of one of the neutrally 
stable linear modes, the nonlinearity causes a wave steepening that leads to the formation of 
weak shock waves and centered expansions on opposite sides of the sheet. The slow 
evolution of disturbances to the vortex sheet occurs on a scale that is large in comparison with 
the wavelength by a factor proportional to the ratio of the wavelength to the amplitude of the 
disturbance. Similar formation of nonlinear "kink" modes for specific incidence angles has 
been observed in numerical calculations by Pedelty and Woodward [3] and Bassett and 
Woodward [4]. 

For a shear layer of nonzero thickness, linear theory leads to an eigenvalue problem 
whose solution gives the complex wave speed. Blumen, Drazin, and Billings [5] have 
shown that an expansion for wavelengths large in comparison with the thickness gives in the 
limit the same neutral wave speeds as for a vortex sheet, subject to the same condition on the 
Mach number. The wave speed, however, has a small imaginary part, implying a slow 
exponential growth of small disturbances in a time or distance larger than the period or 
wavelength by a factor proportional to the ratio of wavelength to shear-layer thickness. 

Balsa [6,7] has observed that nonlinearity influences this shear-layer instability when 
the ratio of disturbance amplitude to shear-layer thickness becomes of the same order as the 
ratio of shear-layer thickness to disturbance wavelength. In a long-wave limit when these 
two ratios are of the same small order, the asymptotic description of slow changes in the flow 
requires matching of solutions for the inviscid external flow, for the main part of the shear 
layer, and for the critical layer where the flow velocity is close to the wave speed. The result 
is a pair of partial differential equations describing a nonequilibrium, nonlinear, viscous 
critical layer, coupled with a condition for the velocity jump across the critical layer in the 
form of an evolution equation for the amplitude. Numerical solutions show the beginning of 
vortex rollup within the critical layer on a spatial scale larger than a wavelength by a factor 
proportional to the ratio of wavelength to thickness. 

In the case studied here, the disturbance amplitude and shear-layer thickness are taken 
to be of the same order. The formulation can be thought of as extending the formulation for a 
vortex sheet by incorporating the effects of nonzero shear-layer thickness. That is, a limiting 
case is chosen such that the thin shear layer experiences slow spatial growth of small long- 
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wave disturbances that are described in terms of a quasi-equilibrium nonlinear critical layer 
and an external flow with nonlinear wave steepening. The problem statement is outlined 
below in Section 2. Solutions for the external flow, for the main part of the shear layer, and 
for the critical layer are given in Sections 3, 4, and 5; a preliminary version of these results 
was given in [8]. Some other limiting cases are discussed in Section 6, corresponding to 
different amplitudes and/or locations along the shear layer. The present formulation and that 
of [6], [7] are seen to represent different special limits of the equations, and a third special 
limit can also be identified. 


2. Problem statement 

A thin laminar shear layer between two parallel streams is subjected to a two- 
dimensional disturbance of small amplitude, traveling in the direction of shear-layer growth, 
with parameters chosen such that the relative disturbance speed is supersonic both above and 
below the shear layer, and such that the growth rate of the disturbance is small. Asymptotic 
solutions are sought in each of three flow regions: the external inviscid flow, the main part of 
the thin shear layer, and the still thinner critical layer where the flow speed and disturbance 
speed are nearly equal. 

There are four independent length scales: a viscous length (based on a reference value 
of flow velocity), a representative shear-layer thickness, a disturbance wavelength (or spatial 
period), and a typical transverse displacement of the shear layer. Three nondimensional 
ratios can then be defined, for example, as follows: £ is the ratio of disturbance amplitude to 
wavelength; 5 is the ratio of shear-layer thickness to disturbance wavelength; and Re is a 
Reynolds number based on the wavelength and on the reference values chosen for the velocity 
and kinematic viscosity, say the values in the undisturbed upper stream. Asymptotic 
solutions for the perturbed flow variables are then sought in the limit as the three parameters 
£, 5, and 1 /j Re approach zero, in terms of suitably defined coordinates corresponding to each 
of the three flow regions. The relative sizes specified for the three small parameters will 
depend on the choice of a particular limiting case to be studied. Here the thickness and 
amplitude are taken to be of the same order, so that 

£ = 0(5) (2.1) 

A relation between the orders of 5 and 1 /Re is chosen in Section 5 so as to provide a balance 
between viscous effects and slow distortion of the flow. 

All variables to be introduced are nondimensional. Lengths are made nondimensional 
with respect to the wavelength (or spatial period) of the disturbance. Velocities and 


3 


thermodynamic variables are made nondimensional with respect to some convenient set of 
reference values, say the values of these quantities in the undisturbed upper stream. 

The coordinates x and y lie along and normal to the undisturbed stream above the 
shear layer, with origin at the origin of the shear layer. Changes occurring on scales 0( 1/e) 
are represented in terms of "slow" variables x, y, and 7 chosen to agree with the scales for 
slow nonlinear wave steepening in the external inviscid supersonic flow: 

x = £(x-L), y = ey, 1 = 3, (2.2) 

where L is a nondimensional downstream distance characterizing the region where nonlinear 
disturbance growth is to be studied; the use of analogous variables in the unsteady one- 
dimensional case is shown in [9]. The relationship between the orders of 8 and 1/Re chosen 
in Section 5 will be such that the slow changes do occur on a scale larger than a wavelength 
but smaller than the length of the shear layer. The slow spatial growth of the shear layer due 
to viscous diffusion will introduce only higher-order changes on the scale considered, since 
the flow is nearly parallel on the scale given by (2.2). Moreover, the spatial extent of the 
assumed disturbance is taken to be small in comparison with the distance from the origin of 
the shear layer, but large in comparison with the scale for the slow variable x\ solutions will 
be sought on the scale of x but not on the still larger scale measuring the extent of the 
disturbance. Temporal or spatial instability can be studied; spatial instability will be 
considered here, since the flow is expected to be convectively unstable [10]. The dependence 
on 7 is retained only through the beginning of Section 3, and is then dropped for spatial 
instability. 

The nondimensional velocity components in the x and y directions are u and v, and c is 
the nondimensional speed of a disturbance traveling in the positive ^-direction. A coordinate 
X is measured relative to this disturbance, and a coordinate Y is measured from a reference 
location y=A within the shear layer: 

X - x-ct, Y = y-A(X£i ) (2.3) 

The nondimensional pressure, density, temperature, and viscosity coefficient are p, p, T, and 
p. The Mach number M is defined in terms of the reference quantities mentioned, as is the 

Prandtl number Pr, which is assumed constant. The nondimensional velocities, temperatures, 
and densities in the undisturbed upper and lower streams are denoted by u=u + , T-T + , and 

p=p + . The undisturbed pressure is taken to be the same above and below the shear layer, 
and it is assumed that the ratio of specific heats is likewise the same; the specific heats are 
taken to be constant. If the reference quantities are the values in the upper stream, then 
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u + =T + =p + = 1 and p + =p_= 1; the notation u + , T + , and p + will, however, be retained for 
convenience. We take u + >u_ and consider disturbances traveling in the direction of 
increasing x, with relative speeds supersonic both above and below. Some of the notation, 
including the expansion (3.5) of A, is indicated in Fig. 1. 


A shear-layer coordinate y is defined by 


y 



(2-4) 


and the undisturbed velocity and temperature profiles are expressed in terms of this coordinate 
by u= u Q (y) and T=TJy). The displacement function A introduced in (2.3) is defined such 
that y = 0 when u Q =c. A critical-layer coordinate y is needed in Section 5, and is defined 


there by 



(2.5) 


The starting point should be the full set of conservation laws, containing all the terms 
representing effects of viscosity and heat conduction. For brevity, many terms are omitted 
here because they are not required in the leading approximations for the description of an 
inviscid flow or a thin shear layer. The needed terms in the continuity, x-momentum, 


y-momentum, energy, and state equations are, at most, 

P t + (P“) x + ( pv) y = 0 (2.6) 

P(u t + UU X + ™ y ) +^P X = ]&(Wy)y + - ( 2 - 7 ) 

p(v, + uv x + vy + -^p y = ... (2.8) 

P <T, + “T x + vT y ) - (p t + up x + vp y ) 

= + fy 2 + - (2.9) 


P = pT (2.10) 

In the external flow the right side of (2.9) is of higher order than the terms needed, and the 
entropy is nearly uniform: 
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p = (p/p+) r + - . 


( 2 . 11 ) 


dp = y (p/p ± ) y 1 d (p/p ± ) + - 


where the upper and lower signs again refer to quantities in the undisturbed upper and lower 
streams respectively. Within the shear layer, dissipative effects can not be ignored, and it is 
convenient to introduce a differential equation for the spanwise vorticity component 

n=v x -u y ( 2 . 12 ) 


Taking the curl of the vector momentum equation, and using the continuity and state 
equations, again including only terms that will be needed for the thin shear layer (as can be 
verified later), we have 


(Q\ 

ra\ 


ra\ 

P - + pu\ 

— 

+ pv 

— 

[pJt 1 

IpJ 

X 

IpJ 


yM 2 p 


V* + h L p w 


yy 


(2.13) 


3. External flow 

In the external flow the length scale in each direction is the wavelength, or spatial 
period. The perturbed flow variables are functions of the coordinates X and Y introduced in 
(2.3), and of the slow variables defined by (2.2), and are assumed to possess expansions of the 
form 


u = u + + eu } (X,Y^x,yJ) + e 1 u 2 (X,Yjc,y,l) + ■■■ 

(3-1) 

v = ev^X.YJy,!) + e l v 2 (X,Y^,y,J) + - 

(3.2) 

P = 1 + £p x {X,Y^,yj) + £ 1 p 2 (X, Yj,y,7) + - 

(3.3) 

— = 1 + £p.(X,Y*,y,J) + f 2 p 2 (X, Y^xjyJ) + - 
P± 

(3.4) 

The displacement function A is expanded as 


A = eAj(X,x,7) + £ 2 A 2 (X,x,7) + ••• 

(3.5) 

All the perturbations are expected to depend on the ratio e/8 as well as on the variables 
shown, and the dependence on 7 will be dropped later. It is in fact premature to specify the 
orders appearing in the second approximation before the solutions for the first approximation 
have been completed in all three flow regions. Once the matching of these solutions has 
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been carried out, each of the functions « 2 > v 2 > p 2 , anc * p 2 will be seen to contain a term 
proportional to ln£. These terms do not play a major role in the solution; for conciseness of 
notation, they will not be shown explicitly. 

The first approximations to the continuity, X-momentum, and y-momentum equations 
are found by retaining only terms of order £, as follows: 


(u ± v j y ~ ® 


(3-6) 


(“± C ^ U \X + p m 2 ^ XX ~ ° 


(u± c)v xx + p M iP\Y ~ 0 


(3-7) 

(3.8) 


Here the quantities u + -c and c-u_ can be thought of as the propagation speeds of a 


disturbance in the shear layer relative to the flow speeds in the upper and lower streams; or, 
u + -c>0 and u_-c<0 are the flow speeds above and below as seen in the X,Y coordinate system 

moving with the disturbance. The two quantities denoted by M 2 p + (u + -c ) 2 are then the 
squares of the relative Mach numbers based on these speeds and on p ± . The parameters are 
chosen such that each of these expressions is greater than one, and thus that the relative 
speeds are supersonic both above and below the sheet. In other words, we consider cases 
where 


u + 


p_ m M 


< c < 


u + - 


p + m M 


(3.9) 


The equation defining c is given in Section 4, and is the same as for a vortex sheet [1]. 


Elimination of u j and Vj above leads to a wave equation for Py 

P±P\XX - P\YY = °’ A± 2 = M 2 p ± (u ± -c) 2 - 1 (3.10) 

where (5 + >0, as noted previously. Solutions for pj, «j, and Vj that describe disturbances 
traveling outward from the shear layer are then found in the form 



P± u \ _ @±P\ 

u ±~ c P+(« + - c ) 2 M 2 


( F l±>£ (£*’'>>) 


where ^ is a coordinate that is constant along linearized outgoing characteristics: 


(3.11) 
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Z = X-p + Y 


(3.12) 


These are the usual results of linear supersonic aerodynamics. While £ has not been written 
with a subscript, it is to be understood that the definition requires (3 + when y>0 and (3_ when 
T<0. The notation F J+ allows the possibility that the function F J+ for Y > 0 might be different 
from the function Fj_ for F<0. It is anticipated, however, that in the first approximation the 

streamline slopes are the same above and below the shear layer as Y— »±0, and that each is 
equal to the leading term in the slope A x of the surface (3.5). Thus, after integration. 


F 1+T=0 ~ F 1-T=0 ~ A \ 


(3.13) 


and therefore also F, and Fj_ are the same functions of the argument £, everywhere, for Y*Q. 
This result will be shown more clearly in Section 4, but for simpler notation will be assumed 
here as well, and the replacement 


(V{ = = 6 < 314 > 

will be made below, where it might again be emphasized that the definition (3.12) for £ is 
different for Y>0 and F<0. 

The terms of order e 2 in the continuity, X-momentum, and T-momentum equations 
give the equations to be solved for the second terms in the expansions of u, v, and p: 

(m_j_— c) P2^- + u 2x V 2Y ~ \xP\y ^ U \P\X + V l^lP 

“ (P\7 + u ±P\x^ + A \X U \Y “ P\( u \x + ~ (“ix + v iy) (3.15) 


(“+" c ) U 2X 


p + M‘ 


2 P2X - ( U ± C ^ix“ir ( U \ U IX + V l M iy^ 


(■ U\ t + u ± u \x> + p M 2 ( ^ 2>) P\P\X + A \xP\Y P\x) 


(3.16) 


(«± ^ m iP2Y ~ ( u ± C ^ix v ir ^ M i v ix + v i v iP 

- (v lT + u ± v lx> + (~ (Y~ 2 ) P]Pi y - Pty) 


(3.17) 


Elimination of u 2 and v 2 leads to an inhomogeneous wave equation for p 2 : 
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fi±PlXX PlYY - (7+1) 


P±( u ±~ cfM 6 
\ 2 




- P+(u+- C ) 2 M 2 (3 ± 2 (2A lx 0^+ A J xx 6£ - 2p ± (u ± - c) 2 M 2 6- - 

~ j-p± 2 (u ± -c?M*6p - jP ± {u ± -c) 2 M 2 (m 2 P± u ± {u ± -c)- l)^ 
p± P+ 

(3.18) 


The perturbations are taken to be periodic in X and bounded as y— >°o. In order that 
the solution for p 2 contain no secular terms that grow indefinitely as £ increases, the terms on 

the right-hand side of (3.18) that depend only on £ (and have no separate dependence on X ) 
should add to zero: 


(Y+l)p+ 2 (u + - cfM 4 I , 

^2 (9e tk - e fy - J^ u ±- c ^ 0 fT 


p+ 


(3.19) 


where d=(F l+ ) ^ , as already noted. This equation describes the slow nonlinear evolution of 

perturbations in the external flow. With the necessary changes in notation taken into 
account, it can be shown that the result (3.19) is equivalent to that of [2] in the special case 
u + =-u_=T + =l. The nonlinear term and the term containing a derivative with respect to y 

are the same as for the far field of the flow past a rigid surface; these terms could have been 
derived instead through matching of an inner solution in terms of y with an outer solution in 
terms of y. The remaining terms allow for slow spatial or temporal distortion of the surface. 
Eq. (3. 19) can be integrated twice with respect to q to give a differential equation for F { . The 
arbitrary function of x and y arising after one integration is zero if Fj is required to be 

periodic in the second such function remains arbitrary, and may be also taken equal to zero, 
since Fj plays the role of a potential function. For consideration of spatial instability, the 

derivative with respect to 7 is omitted throughout the following. 


Equation (3.19), after one integration and with the dependence on t omitted, states 
that 6 is constant along characteristics that can be defined in the form 
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es + x 0 



(3.20) 

where s is measured along the characteristics. On a given characteristic, identified by the 
values of 3c and 0 is equal to its value at the point where that characteristic intersects the 

surface y=0, which in turn equals the value of A lx at that point: 

= <*Vo’ 0) = A ix ( V F o ) (321) 

If A^Xpc) were known for one spatial period, say 0<X<1, and for all 3c>0, where the 

disturbance is regarded as introduced at 3c=0, the differential equation (3.19) would allow 

continuation of the solution away from y=0. For a vortex sheet, the continuity requirement 
for p and the kinematic conditions for v would lead to a differential equation for Aj equivalent 

to a result given in [2]. But in the present case the nonlinear external flow is coupled with 
the flow in a nonlinear critical layer, and so the solutions for the main shear layer and critical 
layer are needed in order to provide the initial values of flat y=0 and x>0. 


A particular solution for p 2 is easily found from (3.18), after the terms shown in (3.19) 
have been omitted; particular solutions for « 2 and v 2 then follow from (3.16) and (3.17), and 
of course will also be consistent with (3.15). Solutions to the associated homogeneous 
equations are chosen to include only outgoing waves, which are represented by undetermined 
functions of £ denoted by F 2+ for 7>0 and F 2 _ for 7<0; it is anticipated that F 2+ *F 2 _. The 

full solutions for p 2 and v 2 , with the dependence on T omitted, are then 

p 2 = -p ± (« ± — c) 2 M\d 4 - (7-1) P ± \u ± - cf M 4 \ fl 2 + -j- p ± (u ± - c) 2 m\f 2± )^ 

P± P± 7 

(3.22) 

V 2 = - P+( u +- c)A , - uJF j— - + (n ± -c)(F 2± )^ (3.23) 

P± 

where Aj=Aj(X^) and 6=F l ^=9{^pc,y); the ± subscripts have now been omitted from Fj. 
The particular solution for p 2 is chosen for later convenience to include the second term, as 
well as the first, in (3.22). The corresponding particular solution for v 2 then contains the first 
three terms shown in (3.23). Expressions for « 2 and p 2 are also easily found. Related results 
for second-order supersonic aerodynamics were given by Van Dyke [11]. 
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If, after (3.19) is integrated twice, F ^ is eliminated from the solution (3.23) for 
v 2 , and /? 2 is found by substituting (3.4), (3.1 1), and (3.22) in (2.1 1), then 

v 2 = -P± (u ±- c)A \ e i, + ~fh F <* 

' + (« ± -cXFa3{ 

P± 

P 2 = -p> ± -c) 2 yM 2 Aje^ + -^-p ± (n ± -c) 2 yM 2 (F 2± )| 

P± 

The term proportional to Fj- agrees, as it should, with the corresponding term given in [7]. 

The remaining terms can be shown to be consistent with the second-order pressure-angle 
relation for isentropic flow turning through a simple wave. For points close to the shear 
layer, i. e., as Y— >0, the expansions for v and p are found to have the form 

v — e{u_£-c)[A^ — P ± YA lxx + ••• ] 

+ c 2 [— (w_j_— c)f3+Ay A + (p+~ c X^2±)x 

+ p 2 A lx ~ 4 p $ P±( u +~ c) 5 M 4 A^ + •••] + ••• (3.26) 

P = 1 + £yM 2 ~-p ± ( u ± -c) 2 [A lx - P+YA }XX + •••] 

P± 

+ E 2 yM 2 ~p ± (u ± -c) 2 [~P ± A l A lxx +(F 2± ) x + -] (3.27) 

P± 

where (F 2± ) x is now evaluated at T=y=0, and 9(Xjc,0) has been replaced by A lx (Xjc). In 

contrast to [2], where a vortex sheet is considered, here the shear-layer thickness is not taken 
to be small in comparison with the displacement of the sheet, and so the expansions (3.26), 
(3.27) are to be matched with solutions obtained for the main part of the shear layer. 


(3.24) 

(3.25) 
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4. Main shear layer 


Within the shear layer the variables u, p, p, and T are assumed to possess expansions 
in terms of the shear-layer variable y=Y/S in the form 


u = u 0 (y) + eu^Xys) + - 

(4.1) 

p = 1 + £Pj(X,y^) + £8p 2 (X,yj) + - 

(4.2) 

P = P 0 (y ) + £p x (X,yj) + - 

(4.3) 

T = T 0 (y) + ef^Xyjc) + - 

(4.4) 


where again the perturbations depend on the parameter £/ 8=0(1) as well as on the variables 
shown. The formulation resembles closely that of [7], which also considers a long-wave 
limit, but certain differences appear because the displacement of the shear layer is now larger 
than in [7], and a displacement term y=A has therefore been subtracted in the definition (2.3) 
for Y. It is convenient to define a transformed velocity V by subtracting certain terms to 
allow cancellation of terms introduced in the differential equations as a result of the 
coordinate shift Y=y-A. Thus we also set 

v - (a - c)A x - £uA- = V = edV 2 (X,yjc) + — (4.5) 

where the continuity equation has been used to show that the expansion of V contains no term 
of order c, this is consistent with the anticipated result (3.13). As in the case of the external 
flow, the orders of subsequent terms are in fact not determined until the solutions for earlier 
terms have been completed. In particular, each of the functions u v p 2 ,p v 7^ and V 2 will be 
found to contain a term proportional to Ine; as in Section 3, these terms will not be shown 
explicitly. 

It is easily seen from the transverse momentum equation that p j is constant across the 

shear layer. Matching with the external-flow solution (3.27) for the pressure then shows that 
the quantity p + (u + - c) /f3 + has the same value, say F, both above and below the shear layer, 

and that 

p x = Y M 2 rA ]X (4.6) 

where F is defined by 

r= = p -<“-- c > 2 (47) 

/> + 
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Combining (4.7) with the definition (3.10) of /? + gives (j0 + -/3_)(l-/? + j3_)=O; the case 
leads to a quartic equation for c: 

- 5 — 5 + - 7 — o = 1 (4.8) 

M 2 p + ( u + -c) 2 M 2 p_(u_-c ) 2 

in agreement with the results cited earlier for waves on a vortex sheet [ 1 ] and for long-wave 
disturbances on a shear layer [5]. As anticipated in (3.9), each term on the left side of (4.8) 
is smaller than one. Explicit solutions for c can be found in special cases [1]. 

The next terms in the expansions of the continuity and y-momentum equations and the 
first approximations to the X-momentum and energy equations are: 


“lX + ^ 2y + y ( u 0 C ^P\X ~ 0 

(4.9) 

(m 0 -c) A xxx + y M 2~hy ~ 0 

(4.10) 

(“0 ~ C ) U \ X + \ U Q + y M 2~$\X ~ 0 

(4.11) 

K~ C ^IX + (“0~ G ~ P\X ~ 0 

7 P 0 

(4.12) 


where « 0 ' and T Q ' denote derivatives with respect to y . The leading terms in the expansion of 
the equation of state are: 

~ ~ p, T, 

PqT 0 = i, P X = z 1 + r 1 (4.13) 

Po T 0 

These equations, along with the energy equation (4.12), have been used to eliminate p j from 
the continuity equation (4.9). 

Elimination of from the continuity and X-momentum equations (4.9) and (4.11) 
now leads to an equation for V 2 that can be integrated to give 



hx 

yM 2 p ± (u ± - c) 2 



P+(“+- cf 

— ? - 1 dy 

P 0 (« 0 - c ) J J 


+ ( G 2±)x 


(4.14) 
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where G 2± =G 2+ (X,x) and is to be determined; the upper and lower signs are to be taken, 
respectively, for y > 0 and ycO. It will be seen in Section 5 that each of the functions G 2± will 
include terms containing a factor ln£. The results for p 2 u j, and Tj are found, from (4.10), 
(4.11), and (4.12), respectively, to be 



(4.15) 



Ml 

yM 2 p ± (u ± - cf 


-P±y + 




f y ( P ± (u ± - cf 

±00 W“ 0 ~ C ? 


\ ) 

1 dy 

J ) 


yM 2 p Q (u 0 - c) 


^2± “o 


(4.16) 


^ 1 =- 


(4.17) 


yM 2 p + (u + - cf 


where P 2± =P 2± (X,x) and is to be determined. The integrals shown here are equivalent to 
those arising in the long-wave derivation of [5]. All the terms in (4.14) through (4.17) appear 
also in the expressions given by Balsa [7]; although the limiting cases considered are 
different, these aspects of the disturbances to the main shear layer are the same. Additional 
terms that appear in [7] are absent here as a result of the coordinate change from y/<5to Y/8. 


Since n 0 '— >0 and f 0 '— >0 as y — »±°°, it is easy to determine the behavior of and fj 

for large y, and then to verify that the solutions for u and T match properly with solutions 
obtained for the external flow. But matching of the solutions for v and p as Y — >0 and y— 
is possible only if certain relationships among the unknown functions are satisfied. With 
the help of (4.5) and (4.6), the shear-layer expansions as y — are found to have the form 


- £(w ± c)A^ x + £^ |^(n ± c)A o A,y + u + A.~ 


P + 


+ e 8 (u^_ c) ( P ± y ) + (G^ + ) + 


(4.18) 
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P = 1 + eyM 1 p ± {u ± -cf~-A lx + £5 (- yM 2 p ± (u ± - C ) 2 A lxx y + P^) + ... 

p+ 

(4.19) 


From comparison of (3.26) and (3.27) for the external flow with (4.18) and (4.19) for the 

8 

shear layer, expressions for -G^iX.x, 0) + A 2 (Xjc) and P 2+ (X,x) are found in terms of 
F^JX, x,0) and A,(X,x): 

“ ( G 2 ±)x + A 2X = C ± A lx + ( F 2 ±)x ~ $± A \ A \XX + 4 ) A \x 

(4.20) 

Pt± = rM 2 £ s r((F^) x - P ± A t A lxx ) (4.21) 


where the constants C + are defined by 

c ± = ^ ( t 2 ~ u±> = _1 + ^ ^ 2 “ !) <4 ' 22) 

and the result /?J3_=1 has been used. Each of the quantities £5 (G 2+ ) x and £ 2 A 2a , represents 
a contribution to the streamline slope in the flow at the edge of the shear layer, and it is seen 
that the two terms appear together in (4.20), suggesting an indeterminacy. As noted again in 
Section 5, the definition of A can be made precise if y=A (i. e., T=0) is identified as the line 
where u-c is exactly zero; then £<5(G 2+ ) represents a term in streamline displacement relative 
to the line y=A. The jump G 2+ -G 2 _, rather than the functions G 2+ and G 2 _ separately, will 
be needed in the derivation of an evolution equation for A j in Section 5. Combining (4.20) 
and (4.21), and using the fact that /? + /?_=l in the final term, 


£ ( G 2+ G 2-^X~( C + C ^ A \x 


S J_ P 2+ p 2- 
£ r yM 2 


(4.23) 


The behavior of the shear-layer solutions for small y is also needed. As y— >0, 
expansion of u Q andr o in Taylor series shows that the integral appearing in V 2 , u v and 7j is 

O(lZy). One way of removing the singular part is shown in [7]; a somewhat different 
procedure is used here, leading to expansions of the perturbation quantities as y — >0 in the 
form 
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(4.24) 

(4.25) 

(4.26) 

(4.27) 


where only those terms that may be needed later have been shown. The subscript c indicates 
that a function of y is to be evaluated at the location of the critical layer, namely at y =0, and 
the integrals J + and K + in u 1 and p 2 are defined by 
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(4-28) 



r p Q (u Q -c) 2 
P±( u ±~ cf 


- 1 


dy 


(4.29) 


The integrands in J + have been modified so as to be integrable at y=0 and also to approach 
zero exponentially as y — (if u 0 ~u + —> 0 and T Q -T + — K) exponentially); the first few terms 
in u ] appear because of these changes. The vorticity (2.12) as y-» 0 is 
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(4.30) 
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An approximate stream function will be useful in Section 5 for obtaining solutions of 
the energy and vorticity equations in the critical layer. In terms of variables X, Y, and x, the 
continuity equation is 

[p(u-c)] x + (fiV) y + £{pu = 0 (4.31) 

For y=0( 1) the solutions derived above show that the third term in (4.31) is 0(£ 3/2 ), i. e., of 
higher order than the first two terms. Thus each of the two mass-flow components p(u-c ) 
and pV can be written in terms of a stream function y/, or y/ = y//8if the appropriate scaling 
for the shear layer is introduced, plus a higher-order correction term: 


p(n-c) = \)f Y + - = + - 

pV = - yr x + - = - S y x + — 


(4.32) 


Substitution of the expansions for p, u, and V into (4.32) allows calculation of derivatives of y/ 
as y— >0: 

Vy ~ P0c u 0c ? + ’ Vx = £ FA IXX + ( 4 - 33 ) 

u 0c 


Thus, as y-^0 

^ = \~Poc“oc'y 2 + °(y 3 >$ 


(4.34) 


5. Critical layer 


5.1. Temperature and vorticity. 

One way of recognizing the scale for the critical layer is to note from (4.34) that the 

~ — 1/9 

first approximation for y/is no longer given by the undisturbed value when y = 0(e ) and 

therefore y/= 0(e). This is the expected scale for the thickness of a nonlinear critical layer, 

as given first by Benney and Bergeron [12]. Thus an appropriate coordinate for the critical 
layer is y=y/e I/2 , as anticipated in (2.5). The orders of the largest flow perturbations in the 
critical layer are determined by matching with the solutions for the main part of the shear 
layer. The asymptotic expansions are then assumed to have the form 

u = c + e^u Qc y + eu x {X,ypc) + ••• (5.1) 

P - 1 + £p x {X,yj) + £dp 2 (X,yj) + - (5.2) 


17 


v = £dV 2 (X,y^) + t? l2 8V 5/2 (X,yJ) + - (5.3) 

T = f Qc + £ m T m (X,y$) + eT x {X,yJ) + - (5.4) 

^ 1/9 A A A A 

P = P 0c + £ p 1/2 (*.;y.*) + £p,(Xj^) + - (5.5) 

1/2 

*2 = « 0 c' + ^-A/2^5* 5 ) + |A<***) + - (5-6) 

^ = ¥ = ¥ 0 (X.ys) + £ ¥ m &’y*) + - (5.7) 


The largest terms in the temperature and the vorticity have been taken as constants, equal to 
the values in the undisturbed flow close to y=0; from the definition (2.12) of the vorticity, the 
largest term in the velocity difference u-c is then linear in y. In the expansion of p, it has 
been anticipated that there will be no term of order r^ 2 , since terms of this order are not 
needed in p for the external flow or the main part of the shear layer. Certain functions 
appearing in the asymptotic expansions for these other parts of the flow will be found to 
contain a factor ln£, since, as seen below, the matching of solutions for the critical layer and 
the main part of the shear layer will require the appearance of a term proportional to ln£ in 
each of the functions G 2± . The factor ln£ will not, however, be present in any of the 

functions shown in the expansions (5.1) through (5.7) above. 


The first two terms in the y-momentum equation (2.8) show that p x and p 2 are 

. A A 

independent of y, l. e., are constant across the critical layer. Matching as y— »±°° and 
y— » 0 then shows that p x is equal to p j, and that p 2 is equal to the value of p 2 at y=0; the 

latter also of course implies that there is no jump in p 2 across the critical layer. From (4.6) 

and (4.25), it follows that 

p x = 7 A/ 2 r A ]X (5.8) 

P 2 + ~ P 2- = r M2rA ixx^ K + - K -^ ( 5 - 9 ) 


Substitution of (5.9) into (4.23) leads to an expression for the jump G 2+ -G 2 _ in terms of 
derivatives of Aj and known quantities. The largest term in the continuity equation (2.6) 

A A 

shows that also V 2 is independent of y, i. e., is constant across the critical layer. Matching 
as y— >±°° and y— K), with the help of (4.24), shows that 


A 



Pqc “ O c 


rA 


\XX 


(5.10) 
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The definition (4.32) of the stream function can be rewritten in terms of \j / , where, as 

before, the stream function satisfies the continuity equation only approximately, because the 
higher-order term £(pu) - in (4.3 1) has not been taken into account. The largest terms give 

= Poc“Oc'y • V'0X = -P0 C '2 < 51 » 

A A ^ 

Thus in the first approximation the streamlines are lines with slopes dy/dX=V 2 /(u Qc 'y). The 
first approximation to the stream function is then 

V'O = 2^0c w 0c^ + ^00 ’ *%) = ~ , r ( A lX _(A lA0 (5-12) 

u 0c 

A A A 

where, for each X, i// Q0 is the value of i// Q at y=0, and (A lx ) m is an integration constant (which 
may be a function of x) defined as equal to the maximum value of A ^ in one period, say for 

A 

0<X<1; this choice is made so that y/ QO <0. Also, it is seen that the derivatives of (5.12) 
match properly with (4.33) as y-»±oo andy— >0. Since A } is periodic in X, the streamlines 

have a cat's-eye structure. The saddle points and centers of the cat's eyes are the points 
where \p QQ , p v and A IX reach their maximum or minimum values respectively; at the 

saddle points and W 00 < 0 elsewhere. For each X, the approximate stream function \pr Q lies in 
the range yr 00 < i/r 0 <oo, with 1/^=0 at the boundaries of the cat’s eyes. 


The first approximations to the energy and vorticity equations (2.9) and (2.13) are 


given by 


A A 


“oc'yV mh + y 2^i/ 2 )0 = 0 


'y 

A 


“Oc'y^m^X + V 2^\n)y y M 2 PlX^m^y - 0 


(5.13) 

(5.14) 


A A 

where, from (5.8) and (5.10), p^ x and V 2 are proportional to the equation of state 
(2.10) and the result (5.13) have been used in obtaining (5.14). Thus the temperature 
perturbation is constant along the approximate streamlines, and the vorticity changes along a 

A A 

streamline as a result of the baroclinic torque. That is, T l/2 and the unforced term in jQ j/2 are 
functions of yr 0 and x, and the solutions have the form 




P 1/2 

P0c 


(5.15) 
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(5.16) 
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1/2 


,2 P \ ( 0 l/2 +)( p + Q\12±(Vq& 


yM 


where the particular solution for Q^ /2 can be verified by direct substitution, with the help of 

A A 

(5.12); the equation of state (2.10) has been used to relate p ]/2 and T^ /2 , and the functions 

A A A A A A 

0,/2±«/^ aHd Q m± (Y 0 s) are still to be determined. The notation 0 ]/2± , Q\/ 2 ± 1S i nten ded 
to emphasize that these are different functions of above and below the cat's eyes; in each 

A A A 

of these regions 0< while within the cat's eyes Wqq- Wq<0- 


The second-order energy and vorticity equations are 


A A 


U 0c' yT lX + V 2 T ly ~ U \( T w)x V 5n( T w)y C m)x 


+ — 'y—Pty + ~r~T — (7\ /2 )aa 

r 0c P 0c Plx *s 2 . RePr p 0c m) yy 


(5.17) 
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(5.18) 


where T 1/2 and Q }/2 are given in terms of 0 1/2+ and Qy 2 + W (5-15) and (5.16); equations 

A A A 

(5.13) and (5.15) have been used in simplifying (5.18); and Pj is related to p x and through 

A A A 

the equation of state (2. 10). Since 0 ]/2± and Q ]/2 ± are functions of i// Q , derivatives of these 
quantities with respect to y are found using dl^y=p^Ti^ c 'y d/dy/^. 

In each of (5.17) and (5. 18) the first two terms on the right-hand side can be rewritten 
and grouped with the terms on the left-hand side. Combining the second terms in the 
definition of the stream function yr with (5.1 1) gives 


U l%X + V 5/2 Voy - u 0c'y W m ) x V 2 ^1/2^ 


(5.19) 
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(This is simply the second approximation to the statement that the stream function is constant 
along a streamline.) Equations (5.17) and (5.18) then can be shown to become 
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+ (5 ' 21) 
pQc U 0c y 

where the notation on the left-hand sides indicates differentiation with y/ Q fixed, and so these 

A 

terms represent derivatives along the approximate streamlines. Since the quantities T, and 

A A yy A A ^ A 

f2j appear only in the combinations ^ \~W\i2^ \I 2 )^iq anc * l ^ e solutions for 

A a A A A A 

T j and will contain terms ¥ 1 / 2 ^ 112 )^ anc ^ ^i/2^i/2^ 0 ’ respectively. Thus the 

A A A A I/O A 

quantities 0 |/2+ and Q l/2± can be thought of as functions of y/=y/ 0 +z ¥\/ 2 + '" rather than 

A 1/9 

y/ Q , and the terms mentioned (with a factor £ ) are then the linear terms in Taylor expansions 

A A A A _ AA 

of & l l2 +(WrX) or Q]/ 2 +(V / > JC ) about y/ =y/ Q . A reformulation of the entire critical-layer 
problem in terms of y/, however, does not seem to offer any substantial advantage, and the 
present formulation in terms of X, y, and x is therefore retained. 


The order of the Reynolds number /?e»l, still unspecified in terms of the small 
parameters e and <5, is chosen at this point such that the diffusion terms are of the same order 
as the other terms appearing in (5.20) and (5.21). Since 5 and £ have been taken to be of the 
same order in (2.1), the choice for Re can be expressed by 

Re- 1 = 0{&) (5.22) 
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and then also {^‘^Re)~ l =0{\) in (5.20) and (5.21). The choices (2.1) and (5.22) specify the 
particular limiting case £— >0, 8-^0, Re—*>° with e/8 and 8^ Re held fixed, for which all the 
physical effects represented in (5.20) and (5.21) are retained at the same level of 
approximation. 

Bearing in mind the expansions (5.4) and (5.6) for T and Q, one can see that the ratio 
of diffusion terms to convection terms in the critical layer is 0[(£ il2 8 1 Re)~ l ], which tends to 
zero in the present limit with 8* Re and e/8 held fixed. The significance of this ratio was 
pointed out by Benney and Bergeron [12] and explored in further detail by Haberman [13]; 
the critical layer is inviscid and nonlinear, or viscous and linear, when the parameter 
A^e 3 l2 8 2 Re)~ l is small or large, respectively. (If the reference length were the shear-layer 
thickness instead of the spatial period, as in [12] and [13], the factor <5 2 would be absent in the 
definition of A.) The present case, for small values of A, is a different limiting case from that 
of [6] and [7], in which A is held fixed in the limit and the amplitude is taken to be e=0(8 2 ). 
The critical layer is thinner in [6] and [7], withy = 0(8) there, and terms representing effects 
of diffusion and slow distortion appear in the first approximation. The result is a pair of 
nonlinear partial differential equations in the fast and slow longitudinal variables and the 
transverse critical-layer coordinate. In the present case the complete problem formulation for 
the largest perturbations requires consideration of the second approximation as well. 

A A 

The additional conditions needed for the determination of 0 1/2 + ^ 2 1 / 2 + (5-15), 

(5.16) take the form of statements about integrals of the right-hand sides of (5.20) and (5.21). 
Since the left-hand sides of these equations are derivatives along streamlines, integrals with 

A A 

respect to X along streamline segments having the same values of T ^ and at both end 

points are zero; the right-hand sides then also have zero integrals along the same paths. 
Outside the cat's-eye streamlines, for y/g>0, the solution is required to be periodic in X. It 

follows that the integrals over one period must be zero. Since lengths have been made 
nondimensional with a spatial period of the disturbance, the integrations may be taken, e. g., 

A 

from X~0 to X=l. For the closed streamlines inside the cat's eyes, it is required that T ^ and 

A 

jQj be single-valued, and therefore the integrals around any closed streamline must be zero. It 
can be seen that the right-hand sides of (5.20) and (5.21) are odd functions of y. This 

A A A A A 

follows because 0 1/2+ and Q 1/2+ are functions of if/ 0 , and if/ 0 is an even function of y ; that 

the first term on the right side of (5.21) is odd follows from (5.20). From symmetry 
considerations, then, an integral over the upper half of a closed streamline is just one -half the 
integral around the entire curve. Thus, not only are the integrals of (5.20) and (5.21) equal to 
zero around a closed streamline, but also are zero separately over the upper and lower halves 
of the streamline. 
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The conditions for y/ Q >0 and y^<0 can be written in a single form if integration limits 
denoted by and are functions of y/ Q having different definitions for ^o>0 and yr 0 <0. 
Outside the cat's eyes, the integrations are to be taken over one period, and so we may take 

W = °, X x (y/ 0 )=1 when y/ Q >0 (5.23) 

For the closed streamlines, X Q and X x are such that 1 ^= 1 ^ at these values of X. The 
definitions are then expressed in the implicit form 

% 0 (X Q s) = $/ Q , y/^X^ = \f/ Q when ^0 < 0 (5.24) 

where we take 0<X 0 <Xj<1. 


Since the integrals of the left-hand sides of (5.20) and (5.21) are necessarily zero, the 
required conditions are expressed by setting the integrals of the right-hand sides equal to zero: 


0 = 


J 


v o 


X, A 

u oc'y 


^ f —^(T m )^dX 

KePr j u -C 1/2 yy 


(?&RePr 


Poc 


u 0c y 


(5.25) 


0 = 


_Jl_ f 

yM 2 J 


q A , 

P ix f 


“oc'y 


A 

r, - 


V\n 




~ ~ i A ^l/2^y 

poc u oc y y 


dX 


Jy 


j 

-J 


“0c'> 




1/2 


‘1/2 0c j. 


0c 


dX + 


1 he f 
^fpRe « J 


Poc 






(5.26) 


where the integrands are understood to be evaluated along streamlines. Each of the equations 
(5.25) and (5.26) contains terms representing effects of diffusion and of slow variation; the 
vorticity equation (5.26) also contains an additional integral associated with the baroclinic 

A A 

torque, which can be rewritten in terms of Ty 2 and p 1 with the help of (5.20). 

A result of this kind was given in another context by Goldstein and Hultgren [14], for 
a stage in the downstream development of disturbances to an incompressible shear layer when 
the appropriate slow variable becomes large. The present formulation is a generalization of 
that given by Goldstein and Hultgren, but in terms of parameter expansions rather than 
coordinate expansions, and with the temperature as an additional dependent variable. 
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For streamlines near the centers of the cat's eyes, the diffusion terms in (5.25) and 
(5.26) can be approximated by substituting (5.12), (5.15), and (5.16), with expanded in 

A A 

Taylor series about the value of X at the center, where i//q Q has its minimum value W Q Q m - ^ 

just the constant and the quadratic term are retained in each case, the integrations over X can 
be carried out, giving differential equations in the variables yr^ and 3c that describe the local 

flow behavior. After this is done, (5.25) becomes 


(®l/2±hc + ( const -) ^(V'o ^00m) (®l/2±^ o j^ o + ® (5.27) 

A 

where the derivatives of 0 J/2+ are evaluated at the center of a cat's eye; a corresponding 
equation is found from (5.26), with the help of (5.20). The results show that one solution for 

A A A 

each of 0 1/2+ and <2 1/2+ has a logarithmic form as y/ Q approaches its value at the center. 
Since derivatives with respect to y/ Q become larger here, it is to be expected that there will be 

a small region in which the diffusion terms should no longer be neglected in comparison with 
the convection terms, and that the solution to an appropriate new set of limiting equations in 
this region could serve to remove the singular behavior. Goldstein and Hultgren [14] show 
this to be the case for incompressible flow and conclude that the logarithmic behavior is not 
sufficient to justify omission of the singular solution. Instead they derive an integral 
condition, a "generalized Prandtl-Batchelor theorem," which shows that the vorticity is not 
constant inside the cat's eyes, and which excludes the singular solution since the integration 
constant multiplying this solution is found to be zero. 


In the present case, supplementary integral conditions for the flow inside the cat's eyes 
are found by following a procedure quite similar to that of [14]. The differential equations 

A A 

(5.17) and (5.18) for T j and can be regarded as vector equations, with 

convection terms in the form (w 0c ' y f) x + (.V 2 f) Y each considered as the curl of a vector 

A A _ A A A 

normal to the streamlines dy/dX = V 2 /(uq c ' y)\ the terms involving u j and V 5/2 are rewritten 
in the desired form with the help of (5.19). The line integrals of these vectors around closed 
streamlines are of course zero, and by Stokes' theorem (or Green's theorem in the plane) the 
area integrals of the remaining terms are zero; these latter terms are the same as those on the 
right-hand sides of (5.20) and (5.21), except for the factor « 0c 'y- For each closed streamline, 

the results can be expressed in terms of line integrals around that streamline: 
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(5.29) 


where, of course, y >0 and y <0 for the upper and lower parts of a streamline, respectively, and 
the first integral in (5.29) can be rewritten with the help of (5.20). 


The conditions (5.28) and (5.29) are analogous to the condition obtained in [14] for 
the vorticity in the incompressible case. With a coordinate change from X,y to X, y/ Q the area 
element is transformed from dXdy to (p Q TiQ c ' y)~ l dXdy/Q, and, as for the vorticity equation 


in [14], the conditions (5.28) and (5.29), except for a constant factor, can be thought of as 
integrals over of equations (5.25) and (5.26). Again as in [14], (5.28) and (5.29) provide 

additional information about the behavior of the solutions near the center of a cat's eye. By 
expanding in Taylor series about the value of X at the center, we can approximate (5.28) 

and (5.29) for streamlines close to the center. As V / Q-*V / QQ m , (5.28) leads to a relation 

A A A 

between the values of (© 1/2 +)j and (®i/ 2 ±^ at the center, indicating that ( 0 1 / 2 ± )j), is 

°a 0 

bounded there, and a corresponding result for G 1/2+ follows from (5.29). Thus it is found 

A 

that the conditions (5.28) and (5.29) require the absence of logarithmic terms in 0 1/2+ and 


A A a 

Q m + as % approaches its value at the center and, moreover, that the derivatives ( 6 >]/ 2 ±)j^ 

AAA 

and (<2 1/2+ )^ remain finite as y/ 0 — > y/QQ m - These conclusions represent the new information 

that is contained in (5.28) and (5.29) but not in (5.25) and (5.26). The situation is analogous 
to that in [14], where, as already noted, the additional condition shows that an integration 
constant multiplying a logarithmic term must be zero. 


Throughout the regions with closed streamlines, the temperature as well as the 
vorticity will not be constant, since, as disturbances move downstream, diffusion effects have 
acted for a limited period of time that is not sufficient to allow the development of an 
equilibrium state for these quantities. The transverse spatial scale for spreading (in 
dimensional form, the square root of the product of kinematic viscosity and time) is the same 

A A 

as the scale for the critical-layer thickness, and so 0 1/2+ and Q 1/2± are still nonuniform in the 
region considered, for x= 0 (l). 
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5.2. Evolution equation. 

In order to obtain an evolution equation for A j, it is necessary to match the 

expressions obtained for the velocity u in the main shear layer and in the critical layer, and 
thus to obtain another expression for the difference G 2 +~ ^2- can be combined with 

A A A 

that given by (4.23). This in turn requires knowing the form of ® 1/2± and Q l/2 ± ^ V'o - * 00 ’ 

As y— >0, the behavior of the solutions for T and Q in the main part of the shear layer is 
found from (4.4), (4.27), and (4.30). If it is assumed that direct matching is possible, the 

A A A 

leading terms in the expansions of the critical-layer solutions for 0 1/2+ and 2 1/2 + ^ y - *±°° 
are found from (5.15) and (5.16) as 
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A A A A i /o A 

where y has been related to y/ Q by (5.12) and we have expanded (y/ Q - Wqq) for large y/ Q . 

This proposed direct matching with solutions in the main shear layer does not contradict 
(5.25) and (5.26); it is assumed here that no additional terms of order i// Q _1/2 or larger are 

needed in (5.30) and (5.3 1), and that the matching can be carried out as proposed. That is, in 
contrast to [14], there appears to be no need for an additional flow region between the main 
shear layer and the critical layer. 


1 /o 

In the critical layer, integration of the terms 0(e IS) in the definition (2.12) for the 
vorticity gives 
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(5.32) 


where the lower integration limit is chosen so that «j=0 when y=0, consistent with the 

definition given earlier for the line y=A as the line where u-c is zero. With the help of (5.30) 
and (5.31), we can now calculate the largest terms in u j as y — »±°°. After addition and 
subtraction of w Qc " y inside the brackets in the integrand of (5.32), the expansion (5.1) for u 

A 

gives, as y— »±°o, 
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and are functions of X andx. The notation is intended to indicate that the finite part should 
be taken; the need for this notation will be seen to disappear in the evolution equation for Ap 

since singular terms will cancel. 


Matching with the expansion of (4.1) asy— X), using the expansion (4.26) for tij and 
setting y=y/£ 1/2 , leads to an expression for G 2± : 
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(5.36) 


Thus, as anticipated, the functions G 2+ , and therefore also the solutions for V 2 , 5j, and 7^ 

given by (4.14), (4.16), and (4.17), will each contain a term with a factor ln£. From (4.20) it 
follows that the functions F 2± will include such a term, and therefore each of the perturbations 

in the external flow will have a term proportional to ln£. On the other hand, (5.32) shows 
that there is no term of this kind in u^\ the other solutions for the critical layer likewise have 
no term containing ln£. 


The result (5.36) can now be combined with (4.23) and (5.9), which were obtained 
from the matching conditions for v and p. The result is an evolution equation for the 
displacement function A } (X^): 
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(c + -cja i5 = -^rV)j 


IX 


- - {(V 7 >ix + ^ + -K)A XX + + (7 + -7_)}* (5.37) 


The integrals 7/ + and 7 + , given by (5.34) and (5.35), are functions of X and x, since they 

“ “ A 

depend on the disturbed temperature and velocity profiles in the critical layer, through 6> 1/2 ± 

A A 

and Qy 2 ±' and on disturbed pressure, through i/z^q and therefore A ^ and p,, according to 
(5.8) and (5.12). The integrals J + and K + , defined in (4.28) and (4.29), depend only on the 
undisturbed velocity and temperature profiles u Q and T 0 =l/p 0 in the main part of the shear 

layer, and are thus independent of both X and x. The terms containing these integrals also 
appear in the evolution equation obtained in [7], but with a somewhat different treatment of 
the singular integrand, as noted earlier. The term on the left-hand side of (5.37) likewise is 
the same as found in [7], whereas the term proportional to A lx arises from nonlinearity in the 

external flow and is absent in [7] because a smaller amplitude is considered there. Each of 
the integrals H + , 7 + , J + , and K + appears in a linear way in the velocity u, and (5.37) is obtained 

by matching solutions for u; thus no products of integrals appear in (5.37). The integrands of 

A A A A 

(H + -HJ and (7 + -7_) involve the combinations ( 0 1/2+ )^+ ( 0 i/ 2 J^ o and Q\/2+ + Q\I2- 

According to (5.30) and (5.31), each of these sums is smaller than Ofy/^ 112 ) as thus 

the integrals ( H + -H _ ) and (7 + -7_) exist. These terms in (5.37) arise from integration of the 
vorticity across the critical layer, and differ from the corresponding terms in [7] because the 
critical-layer solutions are different. 


The evolution equation (5.37) for the displacement function A j is the main result of 
the derivation. Its solution requires an initial condition which specifies A ( for one period in 

_ A 

X at a large negative value of x. The integrals in (5.37) involve the functions 0 ]/2+ and 

A A A 

Q 1/2 +» and so equations (5.25) and (5.26), which contain A j as well as 0 1/2 + ^ Q i/2±’ 
must be solved simultaneously with (5.37). Initial and boundary conditions for 0| /2 + 

A __ ~ 

Q]I 2 + are also needed. It is expected that these functions will approach zero as x — >-<»; the 

~ A 

manner in which they might approach zero is considered briefly below. The forms of 0 1/2 + 

A A 

and Q l/2 + as are given by (5.30) and (5.31). For the closed streamlines inside the 

“ A 

cat’s eyes, the integral conditions (5.28) and (5.29) must be satisfied, i. e., (0 1/? .)fi, and 

A 

(©i / 2 _)a remain bounded as the cat’s-eye centers are approached. As in the analogous 
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incompressible problem, the vorticity, and now also the temperature, are presumed to be 

A A 

continuous at the boundaries fy Q =0 of the cat's eyes, but the derivatives ( 0 J/2 +)ty and 

A ^ 

(0, n )'). might be expected to have jump discontinuities there. For a numerical solution, it 
1/z 

may be necessary, as in [14], to study the thin viscous layers about $r Q =0, following Brown 
and Stewartson [15]; this is not done here. 


Further upstream, in terms of the slow variable x, it is expected that the amplitude is 
smaller, and so Aj is smaller. From the expression (5.12) for the critical layer should be 
thinner, with values of y 2 of the same order as a typical value of Ay From matching with 

A 

the temperature and vorticity in the main part of the boundary layer, the values of @ 1/2± and 

A I/} 

<2 ]/2+ will presumably decrease as Aj . If these suppositions are all correct, then the 
integrals 7 + -7_ and K + -K_ across the main part of the boundary layer and H + -H and 7 + -/_ 

across the critical layer will all remain in the evolution equation (5.37), but the nonlinear term 
proportional to A will be absent in the leading approximation. Since the integrals that 

remain will again lead to slow changes on the scale of the x coordinate, it is necessary that the 
derivative A ^ remain in the equation while A ; decreases. A power-law decay of Aj would 

not permit the needed balance of terms, whereas for an exponential decrease A ] and A j- 

would change at the same rate as 5c-»-oo. Thus it appears that the upstream decay of Aj will 

be exponential in 5c; this is noted again in Section 6. 

The small parameters £, <5, and Re~ ] enter the problem formulation in the combinations 
e/8 and 8* Re, which may be thought of as similarity parameters for this particular limiting 
case. The result (5.37) for the evolution equation, as well as the solutions for the vorticity 
and temperature, may, however, not be uniformly valid for large or small values of these 
similarity parameters. Any such nonuniformity implies the need for a different problem 
formulation in another special limiting case corresponding to a different pair of order relations 
among £, 5, and Re~\ To identify other special limiting cases, it may be convenient to 
replace the shear-layer thickness 8 (nondimensional with the spatial period) in terms of the 
Reynolds number Re and the nondimensional distance L from the origin of the shear layer. If 
the shear layer is regarded as having zero initial thickness at some upstream location, a self- 
similar solution is available, and the thickness is known to grow as the square root of the 
distance. Thus we may set 

8 = L m Re~ m (5.38) 

We might then think of taking ML instead of 8 as one of the basic small parameters. If, for 
example, the viscous length and the wavelength are regarded as specified, the choice of a 
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relation between S and 1 /Re is equivalent to a statement about the distance L. That is, if Re is 
specified, we can think of various possible choices for the position L and the nondimensional 
amplitude £. The results obtained thus far are based on the assumption that b=0{Re ), with 
£ and 8 of order Re~ l/4 . Some suggestions are made in the next section concerning the nature 
of the local solutions when different choices are made for the orders of £ and 8 (or ML). 


6 . Other limiting cases 


6. 1 . Vortex-sheet limit. 

If a disturbance at a location where 8 =0(Re~ ]/4 ) has amplitude £ »<5, it follows from 
setting <5/£=0 in (5.37) that the terms depending on the velocity and temperature distributions 
within the shear layer disappear in the limit. Thus the shear layer appears as a vortex sheet 
in the first approximation, with the evolution equation given by 



(This result can, of course, be found more directly by combining the conditions that the 
pressure be continuous across the vortex sheet and that there be no flow through the sheet.) 
The sign of a in (6.1) depends on the velocity and Mach number in the flow above and below 
the sheet. To illustrate the behavior of the solution we consider a > 0. One special case for 
which a is in fact positive is that in which the temperature is the same above and below, so 
that also p + =p_= 1. For a verification that a> 0 in this case, the three factors appearing in the 

definition of a can be considered separately. First, (4.8) can be rewritten [1] as a quadratic 
equation for {c-(m + +«_)/2} . The solution for c shows that 0 < u + -c < c-u_ and therefore 
(since )3 + /3_= 1) that fi + <l< Thus the relative velocity and Mach number are lower on the 

upper side of the sheet and higher on the lower side. Next, from the second form of the 
definition (4.22) for C + , with the help of the result PJ3j= 1 and the inequalities just noted, it 
can be seen that C + -C_>0 for the present case. Finally, the definition (4.7) for F can be 
rewritten to show that TM 2 =(1+)3 + 2 )/^ + =(1 +jS_ 2 )/)3_. The minimum of (1 occurs at 

{3=1 and is equal to 2, so that r^M 4 is no smaller than 4; since y>l, l-(y+l)r 2 A/ 4 /4 is 
necessarily negative. Thus a >0 if p + =p_= 1; we take a >0 in the following without, of 

course, restriction to only this case. 

If now (6. 1 ) is differentiated with respect to X, it is seen that A remains constant 
along characteristics dX/dx^-txA The solution is given implicitly by 

1 A 
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A lx (X,x) = Mx+ooc A lx (X,F)), WX) = A lx (X, 0) 


( 6 . 2 ) 


where $ is a function of a single argument, and is defined in terms of its form at Jt=0, which 
is considered to be known. Thus (6.2) describes the evolution of the shape of the vortex 
sheet in terms of an initial disturbance represented by Aj^X.O). The solution of (3.19) with 

the proper values aty=0, as well as at x=0, can then be expressed in implicit form by 
evaluating (6.2) at X Q> x Q and combining with (3.21). Substitution for X Q and 3c Q from (3.20) 


allows representation of the argument X Q +a9x Q in terms of 9 and the coordinates x, y. 


Thus 
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Eqs. (6.3) and (6.4), in combination with (3.11) and (3.14), can be used to provide a 
representation of the complete flow field in terms of the initial surface deformation. 


A simple example of a periodic shape A^X.T) is provided by the initial values 

Aj(X,0) = cos X, A lx (X,0) = i>(X) = -sinX (6.5) 


where the spatial period here is taken to be 2% rather than one. Forx>0, (6.2) gives the 
surface slope, and the curvature A j xx (X,T) is found by differentiation as 


Aj XX (X,T) - 


- cos (X + ax Aj^(X,x)) 

1 + oix cos (X + ocx A lx (X,T)) 


( 6 - 6 ) 


This expression becomes infinite when ax =1 and the cosine is equal to -1, i. e., when A 1; ^=0 
and x=a~\ atX=(2n+l)7t, for n=0, ±1, ±2, — . The sheet thus undergoes a distortion, with the 
minima of Aj developing corners such that weak shock waves appear for y>0 and weak 
centered expansions for y<0. The result is consistent with that of [2] for p + = 1, with the 
shocks appearing on the side of the sheet where the relative flow speed u-c is smaller, and the 
expansions appearing on the side where the speed l«-cl is larger. 


The flow field is found from the solution (6.3) for 9, which is proportional to the local 
streamline slope, and which for this case becomes 9= F 1( e = - sin E, where E is again defined 

by (6.4). Above the sheet the flow undergoes compression when the curvature is positive, 


31 


i. e., in the range where cos 2 < 0, while below the sheet compression occurs where cos 2 > 0 
and the curvature is negative. The derivative , proportional to the local streamline 

curvature, becomes infinite at the values of y where its denominator first becomes zero, for 
Icos 2 1=1. The singular behavior arises because the characteristics in the X, Y plane are not 
exactly parallel, an intersection of characteristics occurring at each of the points where F j ^ 

first becomes infinite; a sketch of the wave pattern near the sheet is shown in Fig. 2. Weak 
shock waves form at large distances above and below the sheet as a consequence of the 
coalescing of isentropic compression waves, at the points where characteristics first intersect. 
This is the far-field result for steady flow past a wavy wall, but now the "wavy wall" deforms, 
with increasing or decreasing curvature in portions of each spatial period that contain the 
minimum or maximum values of Aj, respectively. The shocks appear at S=(2n+l)n for 
y>0 and at 2=2nn for y<0, at values of y given by 


(7+1)P + 2 («+-c) 4 M 4 _ (_ M 2 p + u + (u + -c)-l _i 


2/V 


y + a 


x — 


'+ 


COS S, 


(6.7) 


where the values cos 2— + 1 are to be taken for y ^ 0. 


As the downstream distance x increases, the point of formation for each shock wave 
gradually moves downward, toward the sheet for y > 0 and away from the sheet for jkO. 
When ax= 1, shock waves reach the upper surface at the locations X=(2n+l)n , and centered 
expansion waves appear at the same points on the lower surface. With further increase in 
downstream distance, a discontinuity in the slope of the sheet is present at X=(2n+l)7t, and the 
solution for A ]X shows that the magnitudes of these jumps continue to increase until ax=n/2, 

when IAj^l=l immediately ahead of and behind each jump. Thus for \<a x<kI2, and y>0, 
the shock waves at the surface continue to overtake, or be overtaken by, compression waves, 
so that the strengths of the shocks and expansions increase at the sheet: the "kinks" at the 
sheet become stronger. For still larger distances, when ax>n/2, the shocks at the surface are 
weakened by interaction with expansion waves, the flow turning angles at the sheet decrease, 
and the "kinks" become weaker. All disturbances then decay, according to this system of 

equations, as ccx—>°°. Although the effects of changes inside the very thin shear layer are 

ignored in the vortex-sheet limit, it would nonetheless be of interest to consider these higher- 
order effects, taking the function A lx to be prescribed by the external-flow solution for a 

vortex sheet, modified in some way to allow for shock discontinuities; with increasing 
downstream distance, changes within the shear layer may possibly grow large enough that 
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their effect on the outer flow is no longer of higher order. When the value of el 8 is not large, 
the effects of nonzero thickness shown in (5.37) do of course influence the outer flow for till 
x, but the shock-wave behavior is expected to be qualitatively similar to that in the case of a 
vortex sheet. 


6.2. Disturbances at other locations 

The formulation of Sections 3 through 5 is characterized by a relatively thick critical 
layer, with diffusion effects appearing first in a periodicity or singlevaluedness condition. It 
might be anticipated that some other limiting case would correspond to a different balance at a 
location further upstream, where the disturbance amplitude would be somewhat smaller and 
the critical layer somewhat thinner, so that the diffusion terms would appear earlier. In 
particular, a first approximation showing a balance of convection and diffusion terms in 
nonlinear temperature and vorticity equations might be expected. 

One way to start is by considering the behavior of the evolution equation (5.37) as the 
slow variable 3c— >-°o, i. e., by moving upstream from the region discussed in previous 
sections. The amplitude function A ] then decreases, and the quadratic term proportional to 
A lx 2 is expected to become smaller in comparison with the other terms. The slow changes 
in amplitude now are the result of the changes in pressure and streamline slope across the 
main part of the shear layer and the jump in velocity across the critical layer, corresponding to 
the terms 0(8/ e) in (5.37). Since these terms do imply a slow growth, the slow derivative 
should remain in the equation, and so the slow variable originally defined in (2.2) should be 
redefined with scale 5 rather than e. Retaining the slow derivative in (5.37) as 3c— h-°° also 
implies the functional form of A j as x — and requires that A j decrease exponentially. 

However, the anticipated balance of convection and diffusion terms will occur when the 
amplitude is only algebraically small in terms of the parameters, and so the magnitude of x 
need only be logarithmically large, i. e., it is only necessary that -3c = <9(ln Re). The slow 
variable should now be measured from a different reference location, say x=x Q , and redefined 

as 3c = <5(x-jc 0 ). Taking 5(x-x Q )=0(l) and 8(x-L)=0(\n Re), and recalling that L=0(Re* f2 ) 

for the derivation in Section 5, it follows that the relative change in the reference locations is 
given by (Z^-x Q )/L =0(Re~ ]/4 In Re) and is therefore small, so that the relative change in shear- 

layer thickness is also small, and we again have 8^=0(Re~^) as in (5.22). 

Next we can look for a possible asymptotic flow description having downstream 
behavior that matches with the upstream behavior found from the previous formulation. A 
balance between convection and diffusion shows that the critical layer appears when 
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y=0(6< 2/3 ) rather than 0(6 1/2 ) as before. This is the familiar result for a viscous critical layer, 
as given by Lin [16]. Nonlinearity is again indicated through the need for a correction to the 
undisturbed stream function, as suggested by (4.34) evaluated for small y, implying a smaller 
order for the amplitude, namely £ = 0(6* /3 ) rather than 0(6). This is then an example of a 
limit in which the parameter (£ 3/2 6^ 2 /?e) -1 of [12] and [13] is held fixed. Asymptotic 
expansions are introduced, as before, for the external flow, the main part of the shear layer, 
and the critical layer, with various differences appearing because of the smaller amplitude and 
thinner critical layer. The derivations resemble those of Sections 3 through 5 in many 
respects, the greatest difference arising in the now simpler formulation for the critical layer. 
The new representations are presumed to describe a stage in the evolution of a small 
disturbance that precedes the stage considered in Sections 3 through 5. It is suggested that 
this earlier stage exhibits nonlinearity in the critical layer but not in the external flow, and that 
in the later stage the amplitude has grown large enough that nonlinearity enters the description 
of the external flow as well, giving the term proportional to A ^ in (5.37). Since the 

disturbance growth is expected to be exponential, the changes occur within the relatively 
small distance indicated above. 

It seems convenient to express this order-of-magnitude information about the position 
and amplitude in terms of logarithms, defining 

L f =(ln L)/(ln Re), £ f =[ln (l/£)}/(ln Re) (6.8) 

With these definitions, V increases as the downstream distance L increases, but £ 4 decreases 
as the amplitude £ increases. Since L=0(Re ) and, as noted, is close to L, the quantity 
O remains nearly constant and equal to 1/2, even if L is replaced by while in passing from 
the earlier stage to the later stage £^ decreases from 1/3 to 1/4, as can be seen by substituting 
£=;6 4/3 in (2.1) and (5.38). Each of the two stages can then be regarded as represented by a 
point with coordinates V and £*”, namely by the points (1/2, 1/3) and (1/2, 1/4). 

In the formulation given by Balsa [6,7], different order relations were chosen among 
the small parameters. The choice that allowed a balance of convection and diffusion terms 
in the first approximation for the critical layer was, in the present notation, S 5 Re=0( 1) rather 
than 6 4 /?e=0(l). From the interpretation (5.38) of 5as a measure of downstream location L, 
it follows that the present formulation describes flow disturbances at a location further 
upstream and with larger amplitude. In the present case, L=0(Re 1/2 ) and the amplitude is 
0(6), or 00?e“ 1/4 ); in [6] and [7], L=0(/?e 3/5 ) and the amplitude is 0(6 2 ), or 00?e -2/5 ). 
Thus the two representations correspond to different paths of approach to the origin in the £, 
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8, Re 1 parameter space. Or, in the coordinates l) and £* , Balsa's solution [6,7] corresponds 
to the point (3/5 ,2/5) rather than the present (1/2, 1/4) for the later stage. 

We can also consider various other limits, i. e., other choices for the orders of 
magnitude of 8 and £ in terms of Re, corresponding to different points in the V , £* plane. 
Each of the special points (3/5,2/5) and (1/2, 1/4) corresponds to a problem formulation that is 
distinguished in the sense that the proper approximate equations at each of these points is 
"richer" (e. g., [17]) than the equations at any location in a small neighborhood of the point; 
the associated limits might be called "distinguished," in the sense defined, e. g., in [9]. It 
appears that only three points in the plane (for V>0 and £^>0) have this property; the third 
such point is identified below. There are also lines in the plane with the analogous property, 
that along each such line the proper approximate equations are "richer" than the equations at 
neighboring points on either side of the line. Each of these lines will represent a boundary, 
such that some particular effect or effects will be important on one side of the line but not on 
the other. The special points and lines identified are shown in Fig. 3. 

If the amplitude is of the same order as the thickness, then £=0(S) and it follows, with 
the help of (5.38), that 

£ f = (l-L f )/2 (6.9) 

Below this line (f »8) the motion is dominated by nonlinearity in the external flow, and the 
motion is approximated in terms of a vortex sheet, as described in Section 6. 1 . Above this 

line ( e«8) nonlinearity in the external flow does not enter into the leading approximation. 
Clearly the point (1/2, 1/4) of the formulation in Section 5 lies on the line. One end point of 
the line is at (0,1/2), where L=0( 1) and e=0(Re ); here the distance L has decreased to the 
same order as the spatial period, whereas all of the present discussion has assumed that L is 
larger. The other end point is at (1,0), where L=0(/?e) and £=0(1); here the shear-layer 
thickness 8 has increased to the same order as the period, whereas it has been assumed 
throughout that 8 is smaller. 

The limiting forms of the equations are different on different portions of the line (6.9). 
For l)>\/2, and therefore 8 4 Re— >°°, diffusion effects are absent in (5.25) and (5.26), so that 

A _ A _ 

T ]/2 is independent of x and varies with x only because of the term in the baroclinic 

torque; however, a diffusion effect is still present in thin layers bounding the cat's eyes. For 
L^<l/2, i. e., for S 4 Re—>0, the critical-layer thickness is 0(£ 1/2 <5) and diffusion terms become 
as large as the convection terms at the point (311,2/1), where L=0(Re i/1 ) and e=0(Re~ 2n ). 
Since e=0(8), this point corresponds to a limit with 8 1,2 Re held fixed. For still smaller L, 
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i. e., for ZJ<3/7, the energy and vorticity equations continue to show a balance of convection 
and diffusion terms, and the thickness of the viscous critical layer remains 0(<5 I/3 ite _1/3 ), i. e., 
0(SL ~ ]l 3 ). Thus, as L continues to decrease, the critical-layer thickness now increases 

relative to the overall shear-layer thickness, and the velocity jump across the critical layer is 
associated with a smaller term in the critical-layer vorticity. It appears that the second terms 
in both temperature and vorticity are now linear in y, as in the undisturbed motion, and that 

A ^ _ 

the third terms are the unknown functions, matching as jyl— with terms 0(e/y) as y— »0 in 
the main part of the shear layer. The evolution equation remains nonlinear because of the 
quadratic term and also because the integral across the critical layer will involve the 
displacement through the solutions for temperature and vorticity. 

The point (3/7,2/7) is then another special point, with nonlinearity in the external flow 
as well as within a quasi-equilibrium critical layer. The external flow and the main part of 
the shear layer would be described here in essentially the same way as in Sections 3 and 4 
above, whereas the critical-layer equations would be similar to those of [6] and [7], but 
without the slow derivative, and written in terms of a displaced coordinate Y defined by (2.3). 
A detailed study of the equations obtained in this limit, for £=0(8) and 8=0(Re ) might 

be of interest; the associated numerical problem should be considerably easier than that of 
solving the equations in Section 5. 

If instead the disturbance amplitude £ is smaller than the overall shear-layer thickness 
but of the same order as the thickness of the critical layer, then two possibilities arise. A 
nonlinear critical layer, with relative changes in the stream function xjr near y=0 that are not 
small, has thickness 0(e l/2 8), or 0[(£URe) 112 ], whereas a balance of inertia and diffusion 
terms for a viscous critical layer gives a thickness 0[(5 2 /?e) -1/3 5], or 0(L ll6 Re~ m ). 
Equating the orders of the amplitude and the critical-layer thickness then gives, respectively, 
e=0(8 2 )=0(L/Re) and £=0(<5 1/3 /?e _1/3 )=C>(L 1/6 #e~ 1/2 ); the two estimates are the same when 
8=0(Re~ l/5 ), or L=0(Re 2/5 ). An inviscid nonlinear critical layer is found if 8 5 Re-^°°, i. e., 
for L»Re 2/ 5 , and a viscous critical layer if 8 5 Re— >0, i. e., for L«Re 2/5 . In the former 
case, then, 

£ f = 1 — L f (6.10) 

and diffusion effects are expected to be absent from the differential equations describing the 
critical layer, but still to be present in thin layers bounding the cat's eyes. The line (6.10) 
extends from (3/5,2/5) to (1,0). On the other hand, the latter case gives 


and the differential equations for temperature and vorticity contain a balance of convection 
and diffusion terms; as on the line £*=(l-Z^)/2 for 0<ZJ<3/7, the unknown functions appear 
in the third terms in the expansions. The line (6.1 1) extends between (0,1/2) and (3/5,2/5). 
It can be seen that expressions (6.10) and (6.1 1) coincide when L=0(/?e 3/5 ) and £=0(/?e _2/5 ), 
which is the case in [6] and [7]. The two lines (6.10) and (6.11) represent the first 
appearance of nonlinearity, when Re^ ,5 «L«Re and \«L«Re^ 5 respectively: above these 
lines the amplitude is small enough that the disturbances are described by linear theory, but on 
or below these lines the critical layer is nonlinear. 

At least two other line segments also correspond to systems of equations that are 
special in some way. The line e^=(3L^-l)/2, where £/<5->0 such that £= 0[(5 3 ./?e) -1 ], 
connects the points (1/2, 1/4) and (3/5 ,2/5). Along this line the formulation is similar to that 
of Section 5, but the scale for x is 1/8 rather than He. The integrals in (5.25) and (5.26) again 
contain terms representing diffusion and slow variation, but e«8, so that the external flow is 
described by linear equations and the evolution equation (5.37) for Aj therefore no longer 

contains the term proportional to A^ x 2 . The line £^=21^/3 passes through (3/5,2/5) and 
(1/2, 1/3), and intersects £*=(l-L^)/2 at (3/7,2/7). Along the segment between (3/5,2/5) and 
(3/7,2/7) the formulation differs from that of [6] and [7], since (5 2 «£«<5and the critical layer 
is a quasi-equilibrium critical layer: the critical-layer equations show a balance of convection 
and diffusion, but slow derivatives are absent. 

The special line segments divide the plane into a number of regions within which the 
flow development is described by simplified equations obtained by taking suitable 
intermediate limits. (In the case of an incompressible shear layer, Churilov and Shukman 
[18] have identified ranges of parameters for which nonlinearity or viscosity or slow growth 
may be neglected.) In each of the examples considered above, the growth of small 
disturbances with increasing downstream distance is exponential and therefore occurs at a 
nearly constant value of L, i. e., along a vertical line in the e\v plane. In the upper part of 
the plane, as already noted, the disturbances are small enough to be represented in terms of 
the linear theory. As £ increases (and £^ decreases), one or more of the special lines is 
crossed until the line £=5 is reached. 


7. Concluding remarks 

The main purpose here has been to suggest a limiting case that incorporates effects of 
nonzero shear-layer thickness into the description of the weakly nonlinear instability of a 
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supersonic vortex sheet. This limit is such that nonlinearity appears in a critical layer as well 
as in the external flow. Thus the slow changes in amplitude are influenced not only by the 
wave steepening in the external flow, as in [2], but also by the changes in pressure and 
transverse velocity across the main part of the shear layer and the velocity jump across the 
critical layer, as in [6] and [7]. A set of approximate equations is derived for each of the 
three regions, such that asymptotic matching of the solutions should be possible. The 
description of the external flow reduces to that of [2] for the case considered there; effects of 
shear-layer structure are added in the present derivation. The perturbations in the main 
boundary layer are described in essentially the same way as in [7], but nonlinearity in the 
external flow is included here, and the formulation for the critical layer is different from that 
of [7], since terms representing diffusion and slow distortion now appear in the second 
approximation, implying a pair of integral conditions, which are coupled with an evolution 
equation that also differs from that of [7]. 

The equations obtained correspond to specific orders, expressed in terms of the 
Reynolds number based on the spatial period, for the distance from the origin of the shear 
layer and for the amplitude of the small disturbances, and would differ if other choices were 
made. The particular limiting case considered here appears to be just one of three special 
limiting cases. The other two correspond to a nonequilibrium viscous nonlinear critical layer, 
but with linear external flow, investigated in [7], or a quasi-equilibrium viscous nonlinear 
critical layer, with nonlinear external flow, that is mentioned briefly here. Asymptotic 
descriptions of the growth of small disturbances thus are characterized by different sets of 
approximate equations, depending on the location along the shear layer and on the disturbance 
amplitude. 
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Figure 3. Identification of special limiting cases for disturbances at different locations and 
having different amplitudes. 
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